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ABSTRACT 

We developed a magneto-turbulent model for the cosmic ray (CR) electrons 
seen in the radio halo clusters of galaxies. Steady state momentum distribution 
functions of the CR electrons are calculated for given spectra of the turbulent 
Alfven waves. The radio spectrum produced by the obtained CR electron dis- 
tribution is compared to the observed radio spectrum of the Coma radio halo. 
We find that the observed radio spectrum of the Coma cluster is well reproduced 
when the spectral index of the turbulent Alfven waves is ~ 2.8. The obtained 
energy spectrum of the turbulent Alfven waves is steeper than that expected 
from the turbulence theory, suggesting back reaction of the particle acceleration. 
The fitting procedure constraints the amplitude of the turbulent Alfven waves 
as well as the spectral index. Then we estimate the dissipation of the turbulent 
Alfven waves, which is found to be less than the cooling rate by X-ray radiation. 
We suggest that the turbulence which is sufficient for particle acceleration is de- 
veloped in the clusters containing the radio halo. It is most likely that cluster 
mergers create the turbulence and seed relativistic electrons. 

Subject headings: acceleration of particles - galaxies: clusters: general -intergalactic 
medium - radiation mechanisms: non-thermal 

1. INTRODUCTION 

The radio halo in a cluster of galaxies is a diffuse nonthermal (synchrotron) radio emis- 
sion seen in the central region of ~ 1 Mpc in diameter. According to Hanisch (1982), the 
typical luminosity is L t ~ 10 40 — 10 41 erg s _1 (10 MHz < v < 10 GHz), and the spectral 
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index is 1 < a < 2. The typical spatial scale of the radio halos is larger than the optical core 
radius, but smaller than the Abell radius. Giovannini, Tordi, & Feretti (1999) inspected a 
sample of 205 clusters from the X-ray-brightest Abell-type clusters (Ebeling et al. 1996) to 
search for new radio halo and relic candidates. They found only 29 candidates. In addition, 
they found that occurrence of the radio halos and relics is higher in clusters with higher X- 
ray luminosity and higher temperature and confirmed that the positive correlation between 
the absence of a cooling flow and the presence of radio halo. The rarity of the radio halo is 
one characteristic to be considered in its formation scenario. 

In clusters of galaxies with radio halos, nonthermal X-ray radiation due to inverse Comp- 
ton scattering of cosmic microwave background (CMB) photons is expected (Rephaeli 1979). 
Indeed, nonthermal hard X-ray radiation was recently detected in a few rich clusters (e.g., 
Fusco-Femiano et al. 1999; Rephaeli, Gruber, & Blanco 1999; Fusco-Femiano et al. 2000) 
and several galaxy groups (Fukazawa 2001), although their origin is still unidentified. 

Among clusters with radio halo, the Coma cluster is the most informative. The halo size 
is 30' — 40' in diameter corresponding to 600 — 800 kpc. The radio halo region contains 
two bright radio galaxies, NGC4874 and NGC4869. The spatially averaged spectral index 
of the halo, with reduction of contamination by the radio galaxies, is a ~ 1.3. Giovannini 
et al. (1993) presented a spatial distribution of the spectral index between 326 MHz and 
1.4 GHz. They found that the radio spectrum of the central halo region with size of ~ 15' 
is flatter, a ~ 0.8, while in the vicinity of the tailed radio source associated with NGC4869 
is steeper, a ~ 1.8. This means that high energy cosmic rays (CRs) are more abundant 
in the halo region than in the vicinity of the tailed source. The flatness of the radio halo 
spectrum suggests particle acceleration in the intracluster space. This has been pointed out 
by Giovannini et al. (1993). The radio spectrum of the outer region of the radio halo is 
steeper (a ~ 1.8) than that of the central region. 

Kim et al. (1990) studied the magnetic fields in the Coma cluster by using the rotation 
measures (RMs) of background radio sources (QSO and radio galaxies). The observed RMs 
of the background sources seen through the Coma cluster have an excess of ~ 38 rad m~ 2 . By 
adopting a galaxy scale of 10—40 kpc as a correlation length, they found the (electron density 
weighted) amplitude of the random magnetic field to be ~ 1 jiG. On the other hand, the 
comparison between nonthermal hard X-ray emission and synchrotron radio emission gives 
the volume averaged magnetic field strength, if all the hard X-ray is emitted through inverse 
Compton scattering of CMB photons (Rephaeli 1979). Rephaeli, Gruber, & Rothschild 
(1987) obtained a lower limit value ~ 0.11 /xG by using the upper limit on the hard X-ray 
emission from the Coma cluster. The first detection of the hard X-ray emission from the 
Coma cluster by BeppoSAX gives the field strength of 0.15 /iG (Fusco-Femiano et al. 1999). 
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The obtained values are significantly smaller than those derived from RMs. However, it is 
important to note that (1) these estimates assume that the emitting volumes of the hard X- 
ray and radio emission coincide, which might not be guaranteed, and that (2) these estimates 
are affected by the presence of additional infrared background radiation fields (Schlickeiser 
& Rephaeli 1990). These simplification can make the derived field strength smaller than 
that derived from RMs. 

The origin of the CR electrons in the clusters with radio halo is still unclear. Jaffe (1977) 
proposed the primary electron model where the radio halo is an emission by CR electrons 
diffusing away from radio galaxies in a cluster. However, this model contains some difficulties, 
one of which is that the diffusion distance does not seem as large as the radio halo extent. 
Dennison (1980) proposed the secondary electron model where the relativistic electrons are 
produced through decay of charged pions induced by the interaction between relativistic 
protons from radio galaxies and thermal protons in ICM. In this model, however, too much 
gamma-ray emission is produced at least for the Coma cluster (Blasi & Colafrancesco 1999). 

Excess of the high energy CR electrons in the halo region, indicated by the flat spec- 
trum, strongly suggests the particle acceleration in an intracluster space. Models considering 
particle acceleration by intracluster magnetic turbulence were discussed by Roland (1981), 
Schlickeiser, Sievers, & Thiemann (1987), and Petrosian (2001). In this paper, we model 
the radio halo in terms of particle acceleration by the intracluster turbulent magnetic fields. 
The turbulence is assumed to be an ensemble of Alfven waves, and the CR electrons are 
accelerated by pitch angle scattering by the Alfven waves. Distribution functions of the CR 
electrons in energy space are obtained by solving a Fokker-Planck equation for the assumed 
turbulent spectra with various power indices. The calculated radio spectra are compared 
with the observed one to determine the energy spectrum of the turbulent Alfven waves. 

Deiss et al. (1997) proposed a merger shock acceleration model which explains the rarity 
of radio halos. Numerical simulations by Takizawa & Naito (2000) show radio emission is 
luminous while merger shocks accelerate CR electrons but decrease rapidly because of inverse 
Compton cooling after the shocks disappear. The transient feature of this model is consistent 
with the rarity of radio halos. We will discuss the relation between our model and theirs. 

In §2, the basic equations to calculate a radio spectrum are introduced. In §3, the 
calculated spectra are interpreted, and the energy spectra of the turbulent Alfven waves is 
determined so as to reproduce the observed radio spectrum. In §4, the generation mechanism 
of the obtained energy spectra is investigated, and over all turbulence structure is discussed. 
§5 is for summary. 
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2. MODEL 

2.1. Resonant Acceleration by Turbulent Alfven Waves 

Energy of CR electrons is gained through a number of pitch angle scattering by Alfven 
waves (e.g., Wentzel 1974; Skilling 1975; Blandford & Eichler 1987; Berezinsky et al. 1990; 
Terasawa 1991). The resonance scattering condition for an electron of velocity along the 
magnetic field v\\ with a wave of a frequency uja and wave number k is 

c^a = kv\\ + fl e , (1) 

where Q e is the electron cyclotron frequency, and the Alfven velocity relates uja and k by 
uja = kv\. This condition means that the electron travels with the Alfven wave along the 
mean field at a locked phase. The wave number resonant with the CR electron of energy E 
is approximated as 



\k I - -=^- 

"res 

(J>v - V A 

where \x is the cosine of the pitch angle, v ~ c is the CR electron velocity. 
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When \i va/c, \k res \ ~ eB /E, where e and E are the charge and the energy of 
the CR electron, and B is the large-scale-field strength. The energy of the CR electron is 
E ~ m c0 c 2 7, where m e o and 7 are the rest mass and the Lorentz factor of the CR electron. 
For E > 100 MeV, k rcs < 3 x 10~ 12 B^ G cm -1 , where B^ G is the large-scale-field strength 
in units of fiG. In this wave number range, the magnetohydrodynamical description of the 
Alfven wave is valid because k <C k max ~ ty/vi ~ 10 -10 B^q cm -1 , where v\ is the ion sound 
velocity (Holman, Ionson, & Scott 1979). This picture is not valid when \x is much less than 
va/c. When /i ~ 0, the resonant wave number k rcs ~ eB /(m c cyvA) ~ 10~ 8 cm -1 , which is 
larger than the /c max for 100 MeV electrons, and therefore well in the whistler or ion-cyclotron 
range. Steinacker & Miller (1992) discuss the gyroresonance of the electrons with parallel 
transverse cold plasma waves including the Alfven waves, whistlers, and cyclotron waves, 
and derive the Fokker-Planck equation describing the evolution of the electron distribution 
function. Also, Felice & Kulsrud (2001) suggest that scattering at fi ~ is affected by transit- 
time damping effect with long-wavelength waves. These processes should be included in the 
model for the radio halos. For this we have to give energy spectra of all wave modes related 
to these processes to obtain diffusion coefficients. However, it is difficult to assume the 
energy spectra of all wave modes realistically because we do not have sufficient information 
about the generation mechanisms of the waves and the interaction between the waves. The 
aim of this paper is to evaluate the energy spectrum of the Alfven waves accelerating the 
CR electrons producing the radio halo. Then we attempt to do this assuming that pitch 
angle scattering by Alfven waves is dominant process, hoping corrections by other processes 
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near \x ~ is small. However, we must note that if the correction to the present form of 
diffusion coefficient by other wave modes is significant, our evaluation of the energy spectrum 
of the Alfven wave mode would not be valid. In this paper, we adopt a pitch-angle averaged 
diffusion coefficient of the resonant scattering as a reasonably good approximation (Isenberg 
1987; Schlickeiser 1989; Blasi 2000). 

The energy spectrum of the turbulent Alfven waves is assumed to be a single power law 
with an index w > 1: 

P(k) = b^^^(^] \ k <k<k max: (3) 
87r k \k J 

where the applicable range, ko < k < k max , is discussed later; the normalization factor b 
indicates the fractional energy density of the turbulent Alfven waves if k max ^> k Q . In this 
case, the pitch angle-averaged momentum diffusion coefficient is given by 

D pp = (m c0 c) 2 ail w , (4) 

where 

W ~ l ."A, 2-w ,r\ 

w(w + 2) c 

(Isenberg 1987; Schlickeiser 1989; Blasi 2000). Here the relation between 7 and k is 

k = (6) 

and 70 corresponds to k . Also the spatial and momentum diffusion coefficients are related 
by 

pp 

(Melrose 1980; Schlickeiser et al. 1987). 

The transit-time damping effect of fast magnetosonic waves plays a more crucial role for 
lower energy (7 < 1836) electrons (Schlickeiser & Miller 1998; Ragot & Schlickeiser 1998), 
which will emit EUV radiation through inverse Compton scattering of CMB photons and 
possibly contribute in lower energy ends of the radio and hard X-ray emission. Therefore, 
this effect should be considered when broader energy band spectrum models are constructed. 



2.2. Basic Equations 

We simplify the problem as follows. (1) The distribution of the CR electrons is isotropic 
and depends only on the absolute value of the momentum | p \. (2) Spatial averaged dis- 
tribution is concerned. (3) Steadiness is assumed for simplicity; reality of this assumption 
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is discussed later. Thus, the distribution function has a form / = f(p). The effect of the 
spatial diffusion is represented by the escape time formalism in the form of —f(p)/T d , where 
the diffusion time is estimated by 

Up) = |p (8) 

where R is the spatial scale, i.e., the size of a radio halo. 

The Fokker-Planck equation for the present model is, therefore, 

1 d J n 2 d f ■ if 
fd- p \ D ^ P dp- pi - Pf 

-^r + Q(p) =0, (9) 

-Ld 

where pi oss indicates the momentum loss due to synchrotron radiation and inverse Compton 
effect, and Q(p) is the particle injection rate at low energies. For the loss term, we use 

Ploss = - . 4 ? 2 [^mag + U ph ]p 2 , (10) 

om e0 c 

where C/ mag = Bq/(8tt) is the energy density of the large-scale fields, and C/ p h ~ 4 x 10~ 13 
erg cm -3 is the energy density of 2.7 K background photons. For Q(p), the mono-energetic 
injection at 7 = 7; is assumed: 

Q[p( 7 )]=?i%(7)-P(7i)], (11) 
where q± is the injection rate averaged over the radio halo region. 

We rewrite equation (9) in terms of 7 ~ p/(m e0 c), in stead of p, and extract factors 
which are independent of 7 to get 

ct7 z 07 

+4-f{^/}-^P + QM7)] =0 (12) 

where, 

4<tt 

^ = ^zr- r \ u ^ + u ph \. (13) 

orrieoc 

The first term of equation (12) represents simple diffusion in the 7-space, while the second 
term causes the upward advection in the 7-space, i.e., the acceleration, with 'velocity' 7 acC ci = 
(w + 2)a l7 w - 1 . 
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2.3. Characteristic Timescales 

Two other characteristic timescales than T<j are used to normalize the equation: one is 
the acceleration time 

7acccl {W + 2)0,1 

and the other is the momentum loss time 

Ti ~ p^r = — • (15) 

Ipiossl a 2 7 

With the help of equations (4), (5) and (7), the spatial diffusion time (eq. [8]) is rewritten 

as 

Td = ^a l7 - 2 . (16) 

From equations (14) and (16), T<j oc 1/T a . This is because the larger the number of scattering 
is, the higher the acceleration rate is and the shorter the diffusion distance is. Thus, the 
shorter the acceleration time is, the larger diffusion time is. By using these three timescales, 
equation (12) can be rewritten as follows, 

TA ( i 2 \ d 2 f , fTA df 



I I \w + 2] + \tJ' <h 



^-W f) -{i) f+T ' Q =a (17) 

To understand the behavior of the solution, it is convenient to use a critical energy 
m e0 c 2 7 c with which T a = T\. The Lorentz factor for the critical energy is 



7c 



{w + 2)a 1 



a 2 

By using j c , T\/T a and T\/T& can be written as follows 

71 ( 7 x w ~ 3 

y c 



(18) 



T a \7« 



(19) 



T\ _ (w + 2) f va\ 2 f JL\ 1 ~ W _ (20) 



T d 9 \Ra 2 J lc V7c 

The large-scale-field strength and the spatial scale are given for each cluster, and hence a 2 
is regarded as a constant. Therefore, the model free parameters are j c , w, and 

When w — 2, the acceleration time (14) and the diffusion time (16) are independent of 
7. In this case, the Fokker-Planck equation has an analytical solution (Schlickeiser 1984). 
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Schlickeiser et al. (1987) applied an approximated form of this analytical solution to the 
radio spectrum of the Coma halo. For solar flares, analytical solutions of the Fokker-Planck 
equation in general cases with w ^ 2 and finite are discussed in Droge & Schlickeiser 
(1986) and Steinacker, Jaekel, & Schlickeiser (1993). Park & Petrosian (1995) discussed the 
effects of the boundary conditions on the analytical solutions of the Fokker-Planck equation 
in general cases. However, these models cannot be applied to our model in their original 
forms because the synchrotron and the inverse Compton losses are neglected there. 

One of our main purposes in this paper is to search the wider range of parameter space 
than the former work to investigate the radio spectra under various conditions. For this 
purpose, it is convenient to use numerical methods to solve the Fokker-Planck equation 
in general cases with w ^ 2 and the radiation loss. Equation (12) is one dimensional 
diffusion equation and is reduced to the algebraic equation of / by the differencing method. 
The obtained equation is numerically solved by the Gaussian elimination. The boundary 
conditions are df /dp = at p — > and / — > at p — > oo. The comparison between the 
obtained numerical solution and the analytic one for w = 2 suggests that the errors are ~ 1 
% in the present calculation. 



2.4. Radio Flux 

We calculate the radio spectrum from the obtained electron distribution function. The 
volume emissivity is given by 

e„ = J P s (p,u)f(p)p 2 dp, (21) 
where P s is the synchrotron radiation spectrum due to a single electron of momentum p: 

P s(P, v) = ^—T 2 B,x / Ks (77)d77, (22) 



m e0 <r j x 

with x = z//[3eB 7 2 /(47rm e oc)]. The peak of this spectral distribution occurs at a frequency 

v = z/ B 7 2 , (23) 

where z/ B = (0.29 x 3eB )/(A7rm c0 c). This relation (eq. [23]) defines a critical frequency v c 
corresponding to the critical energy m c0 c 2 7 c . 

The radio flux is 
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where / is the source thickness, Source is the solid angle of the radio halo. The radio flux is 
calculated numerically from the obtained distribution function to compare the observation. 
In §3, the radio flux is calculated for the radio halo in the Coma cluster, where we assume 
l(= 2R) = 600 kpc, r2 sourcc = 30' x 30' by adopting the size of the radio halo R = 300 kpc. 

2.5. Model Parameters 

The size of the radio halo R = 300 kpc, the gas number density n = 10~ 3 cm -3 are fixed 
for the Coma cluster, and we examine the case with two field strength B = 0.1 and 1 /zG. 
The model free parameters are v c , w, and g ; . 

The injection energy is assumed to be in the range 100 MeV < Ei < 1 GeV. In typical 
intracluster conditions, the cooling time of ~ 100 MeV electrons is fairly long (~ a few Gyr). 
Thus, it is quite natural that there are abundant electrons with this energy range produced 
by some past active phenomena (Sarazin 1999). The source electrons may be produced 
at shocks in ICM (Takizawa & Naito 2000) or leaked out from radio galaxies, or normal 
galaxies. We choose the value of q\ so that the calculated net flux fits to the observed one. 

3. RADIO SPECTRUM 

3.1. Radio Spectrum and Alfven Wave Energy Spectrum 

Before comparison with the observations, we shall see parameter dependence of the 
solution to clarify the physical properties of the model. Four examples of the calculated 
radio spectrum are taken up in Figure 1 with different indices (a)w = 2, (b)w = 2.8, 
(c)w = 4, (d)w = 6 and with the fixed critical frequency at v c = 10 MHz. The vertical axis 
of the figure is the radio flux in units of Jy, and the horizontal axis is the frequency in MHz. 
We assume B = 1 fiG, and g ; is chosen so that the total luminosity in the 10 MHz— 10 
GHz range becomes L ~ 10 41 erg s _1 for each example. The injection energy is E x ~ 100 
MeV which corresponds to the frequency v ~ 0.05 MHz in the models of w =2 and 2.8, and 
E{ ~ 1 GeV corresponding to v ~ 5 MHz in the models of w =4 and 6. 

The shape of the distribution function is determined through competition between ac- 
celeration and loss. So, the useful parameter is T\/T a (see eq. [19]). The ratio T\/T a is 
proportional to 7 W ~ 3 , and therefore characteristic features of the obtained distribution func- 
tion and, in turn, the shape of the radio spectrum depend on whether w is larger or smaller 
than 3. When w < 3, the injected electrons are scattered up to the critical energy 7 c m c0 c 2 
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because T a < 7] for 7 < 7 C , whereas electrons diffusing out to energies higher than 7 c m c0 c 2 
flow back in the energy space owing to the loss (T a > 71). Therefore, the electrons tend to 
gather around the critical energy / y c m C QC 2 . For w = 2, we see in Figure la that the radio 
spectrum has the peak at v c because of this effect. For w = 2.8, the effect is attenuated by 
momentum space diffusion, but we can still see 'the knee' at u c . 

When w > 3, electrons with energies smaller than / y c m e oC 2 essentially flow toward lower 
energies because T a > 7] for 7 < 7 C . However, a small portion of the electrons is scattered out 
to energy larger than / y c m e oC 2 . Once the energy reaches 7 c m e oC 2 , the electrons are accelerated 
because T a < 7]. Moreover, the higher electron energy is, the more efficiently electrons are 
accelerated. As a result, the spectrum shows a concave shape around ~ v c as seen in Figures 
lc and Id. 

We find that for w > 3 the radio spectrum becomes steeper as w is larger. Figures lc 
and Id show that the spectral index for v > v c is ~ 1 and ~ 2 when w = 4 and w — 6, 
respectively. This is due to strong energy dependence of the acceleration time. The net 
acceleration velocity in momentum space is p ~ p/T a because T a is shorter than 7] for 
v > u c . For a steady state solution, the flow flux in momentum space is nearly constant, 

47rp 2 pf(p) = Anp 3 ^- ~ constant. (25) 

Hence, from equations (14) and (25), we obtain 

/(p)«P" (w+1) . (26) 

By using the monochromatic approximation, when the electron distribution function has 
power law index T = w + 1, the radio spectrum follows 

S„ oc z/-( r - 3 )/ 2 = (27) 

This relation explains the radio spectrum becomes steeper as w is larger, and S v oc v~ x and 
oc v~ 2 when w = 4 and 6, respectively. 

When w = 3, T a and T\ have the same 7-dependence, and therefore the ratio T\/T a is 
constant at any 7 so that there is no critical frequency u c (i.e., 7 C ): we find that f(p) and 
the radio spectrum follows a single power law. The larger the value T\/T a , the flatter the 
spectral index is. When T\/T a ~ 1, the spectral index is a ~ 1 in the 30 MHz— 5 GHz range. 
When w — 3, the static solution to the Fokker-Planck equation without the diffusion loss 
term and the source term is / oc p~ a ' 2 / ai = p- 57 V T i (Eilek & Henriksen 1984). For T a /T\ ~ 1, 
/ oc p~ 5 , and in turn, S u oc u^ 1 by using the monochromatic approximation. 
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If w 7^ 3, the critical frequency v c is observable in a radio spectrum and characterize a 
radio spectrum. Determination of v c and the spectral index allows us to evaluate the energy 
spectrum of the turbulent Alfven waves as shown in the next subsection and §4. 

Table 1 shows the result of the parameter survey in the v c — w space, where w is in 
between 2 and 6, u c between 0.01 MHz and 10 GHz, B = 1 p G, and E\ = 100 MeV. When 
a falls in the observed range 1 < a < 2 in between 30 MHz and 5 GHz, the values are given 
in the columns. The model spectra which is too steep with a > 2 or which is too flat with 
a < 1 are represented by 'S' or 'F', respectively. When the spectrum has a peak in the 
observed frequency range, the model is denoted by 'P'. When w > 4, the lower bound of the 
critical frequency is not known. 

We also calculate for the models with B = 0.1 fiG. The relation between a and the 
model parameters are not altered significantly. 

3.2. The Energy Density of Turbulent Alfven Waves 

The value bBl/{8it) means the energy density of the turbulent Alfven waves for accel- 
eration of the CR electrons radiating at v < 5 GHz, when we assume ko = /csghz i n equation 
(3). The wave number A; 5GHz ~ 9 x 10~ 15 B 3 J* cm 1 is defined as the wave number resonat- 
ing to the CR electrons radiating at v — 5 GHz, which is similar to the observed maximum 
frequency. Figures 2a and 2b show the relation between b and w with B = 1 /zG and 0.1 
fiG, respectively. The solid line and the short-dashed line indicates, respectively, the lower 
and upper bound of b, which corresponds to the range 1 < a < 2. 

The upper limit of the energy density of the turbulent Alfven waves is the energy density 
of the large-scale-field Bq/(8h). In figures 2a and 2b, the long-dashed line represents the 
upper limit 6=1. The shaded region indicates the allowed range of b and w. When w < 4 
with B = 1 fiG, b < 1, i.e., the condition is met. However, when w > 4, the range of the 
parameters is constrained by the upper limit. When B = 0.1 /iG, the condition is met in 
the models with w < 3.5. In the models with b > 1, the energy density of the turbulent 
Alfven waves exceeds Bq/(8tt). Such models are unrealistic. 

As noted in the previous subsection, when w > 4, the lower bound of the critical 
frequency is not known. The values of the lower bound of the critical frequency determined 
by the constraint b = 1 are 3 x 10 -6 MHz and 0.1 MHz for w = 4 and 5, respectively, when 
B = 1 fiG. When u c is less than the lower bound, the energy density of the turbulent Alfven 

waves exceeds Bq/(8tv). 



3.3. The Effect of Spatial Diffusion 

The ratio T\/T& (eq. [20]) is rewritten as follows, 

^~3xlO- 2 ( W + 2)( I ^j , (28) 

where we use B = 1 //G, R = 300 kpc, the gas number density n = 10~ 3 cm -3 , and 7 C = 10 3 . 
Also, the ratio T a /T d is written as 




The spatial diffusion term does not play a significant role in most cases because T d is much 
longer than T a and T\ within the observed frequency range (7 > 10 3 ). However, this is 
not the case when w < 2, which may be included in our parameter survey, e.g., w = 5/3 
for Kolmogorov. If, for instance, w — 2, then Xi/T d ~ O.l^/lO 3 ) -1 and (T a /T d ) ~ 0.1. 
Therefore, T d is comparable to 7] at a energy range in which we are concern, and to T a at 
all energies. 

In §3.4 we search for the best fit model for the Coma radio halo and obtain the best 
model at w — 2.8, v c = 10 MHz with B = 1 /iG. The spatial diffusion time of this best 
model is T d ~ 10 11 yr at 7 = 6.4 x 10 4 corresponding to 5 GHz. Since the spatial diffusion 
time is longer than the Hubble time, the distribution function does not reach to the steady 
state. However, the acceleration time and the energy loss time are shorter than the Hubble 
time, and therefore the distribution function becomes similar to the steady state solution 
quickly but only its normalization is different. The obtained physical quantities other than 
the injection rate are fairly valid. The injection rate q\ would be underestimated. 



3.4. Comparison to the Observed Spectrum 

We search for the best fit model for the Coma radio halo by the method of least squares. 
The average value of squared residuals is calculated as 

J = -!- ]T(log 5 cal - log S ohs )\ (30) 

1 'obs 

where S c& \ is the calculated radio flux and S \, s is the observed one drawn from Deiss et al. 
(1997) and the number of the data points A^ obs is 12. The radio flux (eq.[24]) is calculated 
with the solid angle Source = 30' x 30' and the source thickness I (= 2R) = 600 kpc. (The 
distance to the Coma cluster is assumed to be 70 Mpc by adopting the Hubble constant 
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H = 100 km s" 1 Mpc -1 .) The obtained values of J with B — 1 fxG are presented in Table 
2. The best fit model (J = 0.018) is obtained when v c = 10 MHz and w=2.8. Figure 3a 
shows the radio spectrum of this model with B = 1 /zG, and some physical quantities are 
summarized in Table 3. For the best fit model, b ~ 1.7 x 10~ 5 and the energy density of the 
turbulent Alfven waves for the acceleration is ~ 7 x 10~ 19 erg cm~ 3 . (The factor b represents 
the fractional energy density of Alfven waves accelerating the CR electrons radiating at v < 5 
GHz [cf. §3.2].) This required energy density is much less than the energy density of the 
large-scale-field in the cluster Bq/(8it) ~ 10~ 14 erg cm~ 3 . 

Also, we obtain small values of J when w = 4.5 and v c < 1 MHz. The radio spectrum 
with w = 4.5 and v c = 1.0 MHz is shown in Figure 3b. The required injection rate q (= 
Anpfqi) ~ 3 x 10~ 9 cm~ 3 s _1 becomes very large when E\ = 100 MeV. An injection at higher 
energy might resolve this difficulty. 

Figure 3 shows that the spectra with a small value of J ~ 0.02, can have different 
shapes. If the spectrum shows exponential cut off such as shown in Figure 3a, the index w 
is smaller than 3, or otherwise, if the radio spectrum has a concave shape, the index w is 
larger than 3. When w = 3, the radio spectrum follows a single power law as mentioned in 
§3.1. Deiss et al. (1997) suggests that the value of the spatially integrated flux above 1.4 
GHz is underestimated because of the small integration area. If more diffuse radio power is 
observed above 1.4 GHz, the model with w > 3 may well explain the observation. Further 
observation to determine the radio spectrum at higher frequencies set a stronger constraint 
on the energy spectrum of the turbulent Alfven waves. 

When B = 0.1 /iG, the models with (w, z/ c [MHz]) = (2.8, 10), and (4.5, 1) reproduce 
the observed spectrum. Table 3 shows that the required fractional magnetic energy density 
b is larger than that in the case of B = 1 /zG. When w = 4.5 and v c = 1.0 MHz, the energy 
density of the Alfven waves for acceleration exceeds -Bq/(87t) as noted in §3.2. 



Let us consider physical status of the turbulent Alfven waves in our best fit model. 
Ruzmaikin & Shukurov (1982) gives the rate of energy-transfer via the nonlinear interaction 
among different wave numbers as 



4. 



DISCUSSION 



4.1. The Energy Spectrum of Turbulent Alfven Waves 




(31) 
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where p is the mean gas density. If the energy spectrum of the turbulent Alfven waves is 
determined only through this energy cascade process, it becomes P(k) oc k~ 3 l 2 in a steady 
state because (dE/dt) n \ is constant for k (Ruzmaikin & Shukurov 1982). However, this is 
significantly flatter than that in our best fit model whose exponent is —2.8. Thus, other 
physical processes might work to steepen the Alfven wave spectrum. Obviously, the Alfven 
waves lose their energy through the particle acceleration. The CR electrons between p and 
p + dp absorb the energy of Alfven waves at the rate, 



where we use dp ~ p, 7 aC cei = l/T a . We checked whether this plays a significant role or not in 
our best fit model. For simplicity, we approximate P{k) oc k~ 3 and f(p) oc p~ 6 . In this case, 
(dE/dt) n \ oc k~ 3 and (dE/dt)^^ oc k~ x because the resonant wave number is proportional 
to p~ x . Therefore, the energy loss due to the particle acceleration is more significant in 
higher wave numbers and may be attributed to the spectral steepening. We checked if both 
terms are comparable to each other in the observed frequency range (30 MHz — 5 GHz), 
which corresponds to the wave number range of (10- 14 ~ 10- 13 )5j? cm" 1 (eq.[6]). The 
critical wave numbers k t where (dE/dt)^^ = (dE/dt) n \ are shown in Table 3 for our best 
fit models. They are certainly within or near the wave number range mentioned above. The 
obtained steep wave spectra, w ~ 2.8 and 4.5 for the Coma cluster and w > 2.0 in general, 
could be brought about by this effect. We suggest that time dependent and self-consistent 
calculations for non-linear wave-particle systems with back reaction of particle acceleration 
are necessary to investigate evolution of the systems. 

Next, we make a crude estimate of P(k) in k < k t and check the consistency of our 
consideration mentioned above. We adopt P(k) of our best fit models for k > k t and assume 
P{k) oc k~^ for k < k t . As a result, the broader band model power spectrum is, 



where P t = 5 t (-Bo/87r)(u> — l)/k t and k m i n is the wave number corresponding to the scale of 
the cluster, k min ~ 10~ 24 cm -1 . Then, the total energy density of Alfven waves M tot is given 




acccl 



Aixp 2 (m c0 c 2 )7 acccl / (p) dp 

4n(m c0 c 2 )p 3 ^f(p) 
47rm 4 c 5 ( W + 2)a l7 tu+2 /b), 



(32) 




(33) 



by 




w — 1 \ k 



w-€ { k 




(34) 
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where we assumed £ > 1. If we assume M tot is as large as the energy density of the large- 
scale-field Bq/(8tv), then we obtain £ with given k t and b t . When B = 1 /iG, fc t ~ 10~ 14 
cm" 1 and b t ~ 10~ 6 , we have £ ~ 1.6, which is close to the index of the energy spectrum, 3/2, 
formed by the non-linear interaction. Therefore, this result is consistent with the previous 
picture of the non-linear interaction for k < k t and the acceleration reaction for k > k t . 
When B = 0.1 /iG, k t ~ 10" 15 cm" 1 and b t ~ 10~ 3 , we have £ ~ 1.2. 

The rate of the energy transfer from /c3omh z to higher wave numbers is estimated as 
e ~ [-P(^30mhz)] 2 ^3omhz/(/ w a)- Table 3 shows that the values of e are less than the thermal 
X-ray cooling rate except for the unrealistic case with Bq = 0.1 /iG, w = 4.5 and v c = 10MHz. 
Therefore, even if the energy of the turbulent Alfven waves goes into the thermal energy of 
the hot gas by the dissipation at higher wave numbers, the hot gas could not be heated up. 



4.2. The Origin of the Turbulence 

The energy transfer rate Eq at k m i n ~ 10~ 24 cm -1 is estimated as Eq ~ Pq k^ in / pva ~ 
10~ 31 -B^ G [erg cm -3 s -1 ], where Pq ~ -Bg/(87rA; m i n ). Thus, an energy input rate higher than 
Eq is necessary to excite and maintain the turbulence. It is certain that motion of galaxies 
through the intracluster medium can excite turbulent eddies. In this case, however, the 
energy input rate of the steady turbulence is 

einput ~ pvl/k ~ 10~ 33 erg s _1 cm" 3 , 

where l ~ 10 kpc is a dimension of the galaxy, and v ~ 3.2 km s -1 is the turbulent velocity 
at l (Goldman & Rephaeli 1991). This input rate is insufficient for the required e . If 
the energy of the turbulence is input by the system of scale R ~ 1 Mpc with the velocity 
comparable to the sound velocity v s ~ 1000 km s -1 , the energy input rate is 




This input rate is sufficient for the radio halo formation. A merging event between sub- 
structures can input the kinetic energy into the turbulence at this high rate. Moreover, 
cluster merger can produce seed CR electrons through 1st order Fermi acceleration at shocks 
in ICM (Takizawa & Naito 2000). Indeed, recent high resolution MHD simulations of cluster 
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mergers (Roettiger et al. 1999) show that the bulk flow is replaced by turbulent motion in 
the later stages of the merger. 

4.3. The Supply of Cosmic Ray Electrons 

We calculate the total number of the CR electrons N C r from the luminosity and find 
N C r ~ 10 60 for the Coma radio halo, where we adopt the luminosity L t ~ 10 41 erg s _1 (10 
MHz < v < 10 GHz), the spectral index a ~ 1.3, and the magnetic field strength B ~ 1 fj,G. 

We examine how much the radio galaxy can supply the CR electrons. NGC4874 is one 
of the two dominant galaxies in the Coma cluster. The observed flux of the radio source 
associated with NGC4874 is 5 408 mhz ~4x 10~ 24 erg cm~ 2 s" 1 Hz" 1 at 408 MHz (Jaffe & 
Perola 1974). From the observed flux we calculate the number density of the CR electrons 
at the injection energy ~ 100 MeV. For simplicity, we assume that the shape of the source 
is a sphere with a radius r ~ 10 kpc, and that the radio spectrum follows a single power law 
with an index a ~ 0.5. We also assume the magnetic field strength is 10 fiG. We find the 
number density ucr ~ 6 x 10~ n cm ~ 3 at 100 MeV. The injection rate of the CR electrons 
leaked from the radio source is given by ncR ~ 47rr 2 ncRfD, where vb is the diffusion velocity 
of the CR electrons. Taking v d ~ va ~ 700 km s _1 gives an upper limit of the injection 
rate n CR < 4 x 10 43 s -1 . When n CR ~ 4 x 10 43 s -1 , the CR electrons in the radio halo are 
supplied in iV CR /n CR ~7x 10 8 yr. 

The merger shock acceleration supplies the CR electrons in the radio halo region. Tak- 
izawa & Naito (2000) assumed that the total kinetic energy of accelerated electrons is 5 % 
of the viscous energy, which is nearly equal to the energy of the shock heating, and they 
found that the luminosity of the synchrotron radio emission (10 MHz— 10 GHz) is ~ 10 43 
erg s -1 . We find the number of the CR electrons calculated from the luminosity is 5 x 10 62 
and sufficient for the observed radio halo, where we assume the magnetic field strength is 
0.1 fiG and the radio spectral index is 0.7. According to Takizawa & Naito (2000), the radio 
luminosity decreases after the most contracting epoch. We suggest that the turbulent Alfven 
waves reaccelerate the CR electrons injected by the merger shock and maintain the radio 
halo. 



5. SUMMARY 

We have examined whether the particle acceleration by the turbulent Alfven waves 
produces the radio halo in the cluster of galaxies. As a result, the energy spectrum of the 
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turbulent Alfven waves reproducing the Coma radio spectrum has w ~ 2.8. It is most likely 
that the turbulence with sufficient energy for the radio halo formation is originated from some 
merging events between sub-structures. In our steady model, evolution of the turbulence and 
the magnetic fields with the cluster evolution is not considered. By considering effects of 
the cluster evolution, difference between the halo clusters and the non-halo clusters would 
be clarified. Such difference may also be indicated by observation of intracluster hard X-ray 
and thermal X-ray. 

We thank T. Terasawa, M. Tosa, and M. Hattori for helpful suggestions. We would also 
like to thank the anonymous referee for his/her very useful comments. 
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Fig. 1. — Calculated radio spectrum with B = 1 /iG, R = 300 kpc, v c = 10 MHz and 
various values of w. (a)w = 2; (b)w = 2.8; (c)w = 4; (d)w = 6. The radio spectrum shows 
convex and concave shape for w < 3 and w > 3, respectively. 



- 21 - 




Fig. 2. — Relation between b and w. The solid line and the short-dashed line indicates, 
respectively, the lower and upper bound of b, which corresponds to the range 1 < a < 2. 
The long-dashed line represents the upper limit b = 1. The shaded region indicates the 
allowed range of b and w. In the models with b > 1, the energy density of the turbulent 
Alfven waves exceeds Bq/(8tv). 
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Fig. 3. — Two examples reproducing the observed radio spectrum, (a) B — 1 fj,G, v c = 10 
MHz, w = 2.8; (b) B = 1 fiG, v c = \ MHz, w = 4.5; The observed intensity is drawn from 
Deiss et al. (1997). 
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Table 1. The calculated spectral index with w and u c . 



V c w 



(MHz) 


2.0 


2.6 


2.8 


3.2 


3.5 


4.0 


5.0 


6.0 


10 4 


P 


P 


F 


S 


S 


S 


S 


S 


10 3 


P 


F 


F 


2.07 


S 


S 


S 


S 


10 2 


1.02 


0.97 


1.03 


1.35 


2.61 


s 


S 


S 


10 


S 


1.90 


1.44 


0.96 


1.16 


1.79 


S 


s 


1 


s 


S 


2.03 


F 


0.88 


1.08 


1.54 


s 


0.1 


s 


S 


S 


F 


F 


1.02 


1.50 


2.00 


0.01 


s 


s 


S 


F 


F 


1.00 


1.50 


2.00 



Note. — Values in columns are a between 30 MHz and 5 GHz. 
Spectra of a > 2 and spectra of a < 1 are represented by the charac- 
ter 'S' and 'F', respectively. Spectra which has a peak in the observed 
frequency range is represented by 'P'. 
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Table 2. The calculated values of J with w and v c . 



(MHz) 2.0 2.6 2.1 



3.2 



3.5 



4.0 



4.5 



5.0 6.0 



10 3 


0.71 


0.49 


0.30 


0.33 












10 2 


0.14 


0.11 


0.075 


0.053 


0.91 










10 


1.4 


0.14 


0.018 


0.11 


0.070 


0.17 








1 




1.6 


0.20 


0.18 


0.13 


0.063 


0.027 


0.036 




0.1 






1.0 


0.22 


0.16 


0.077 


0.028 


0.030 


0.20 


0.01 








0.24 


0.18 


0.082 


0.028 


0.029 


0.20 
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Table 3. Model parameters reproducing the observed radio spectrum. 







bB%/(8n) b 


(cm 3 


B-l) 


kt d 










(w, f c [MHz]) 


b a 


(erg cm -3 ) 


B i= 100 MeV 


E i= l GcV 


(cm- 1 ) 




h c 




(erg cm -3 s -1 ) 


B = 1 fJ-G 




















(2.8, 10) 


1.7 x 10~ 5 


6.8 X 10~ 19 


3.8 x 10~ 29 


2.3 x 10~ 30 


5.7 x 10" 


14 


6.1 x 10" 


-7 


2.2 x 10~ 33 


(4.5, 1.0) 


1.5 x 10~ 2 


6.0 x 10~ 16 


2.9 x 10~ 9 


4.1 x 10~ 29 


1.5 x 10" 


13 


8.1 x 10- 


-7 


1.8 x 10~ 30 


B = 0.1 




















(2.8,10) 


1.6 x 10- 1 


6.2 x 10- 17 


1.5 x 10- 30 


1.4 x 10- 30 


4.6 x 10- 


15 


1.0 x 10" 


-3 


5.7 x 10- 30 


(4.5,1.0) 


140 


5.5 x 10- 14 


1.1 x 10™ 


2.5 x 10~ 26 


8.9 x 10" 


IS 


8.3 x 10" 


-4 


5.0 x 10~ 27 



a The value b is the fractional energy density of Alfven waves above fe 5GHz = 9 X 10- 15 B^ 2 cm . 
b The energy density of Alfven waves above fcsGHz- 

c The values of the injection rate q(= 4np?qi) for the steady state solution to reproduce the radio intensity. 

d The value kt is the wave number at which the energy loss rate due to the acceleration overcomes the energy transfer rate via 
the nonlinear interaction. 

e The value bt is the fractional energy density above kt- 

f The value e is the energy transfer rate at A:30MHz- 



